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I.  REVIEW  OF  TASK  OBJECTIVE  AND  ESSENTIAL  TASK  ELEMENTS 


OBJECTIVE 


The  proposed  two-phased  effort  is  intended  to  result  in  an  analytical 
procedure  which  is  capable  of  predicting  the  response  associated  with  a 
fender/vessel  interaction.  As  part  of  Phase  I  efforts,  performance  algorithms 
would  be  determined  relating  energy  and  load/deflection  characteristics  for 
various  generic  types  of  fendering  systems.  These  systems  are  typically 
rubber,  pneumatic,  foam-filled,  etc.  and  have  been  investigated  herein.  The 
resulting  algorithms  would  be  used  to  characterize  a  particular  generic  type 
fender  which  would  be  represented  in  the  response  vessel/fender  interaction 
procedure  developed  in  Phase  II. 

The  response  analysis  to  be  developed  would  be  capable  of  estimating  the 
following: 

a.  The  maximum  energy  absorbed  by  a  generic  type  fender  represented  .in 
the  response  problem. 

b.  The  maximum  reaction  load  input  to  the  pier  and  vessel  hull  during 
response. 

c.  The  maximum  local  deflections  occurring  in  the  fender  and  hull  during 
response. 

d.  The  relative  amounts  of  energy  stored  in  fender,  local  hull,  hull 
mode  motions  and  hydrodynamic  dissipation  during  response,  thus 
defining  the  energy  storage  requirements  for  the  fender  system. 

The  benefits  of  fender/vessel  response  capability  would  include: 

a.  The  ability  to  simply  explore  a  specific  vessel/fender  response 
subject  to  design  constraints  such  as  maximum  permissible  reaction 
loads  and  deflections  or  fender  system  energy  absorption  requirements. 

b.  The  ability  to  plug  in  and  out  alternative  fender  performance  charac¬ 
teristics  via  algorithms  and  explore  overall  vessel/fender  response. 

c.  The  capability  of  optimizing  a  specific  vessel/fender  response  for 
a  given  set  of  problem  constraints. 
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ESSENTIAL  TASK  ELEMENTS  -  PHASE  I 


Task  1 

Task  2 

Task  3 


Task  1 


Task  2 
Task  3 
Task  4 
Task  5 
Task  6 
Task  7 


-  Acquire,  review  and  assess  the  performance  characteristics  of  currently 
used  fendering  systems  in  the  available  literature  in  order  to  estab¬ 
lish  a  data  base  of  energy  absorption  and  load  deflection  data. 

-  Determine  algorithms  which  quantify  fender  system  performance  for  generic 
fenders.  Identify  and  rank  energy  absorption  mechanisms  for  these  fender 
systems. 

-  Based  on  the  literature  search  performed,  identify  an  approach  leading 
to  the  development  of  a  rigorous  analytical  technique  for  predicting 
vessel/fender  interaction.  This  technique  will  be  fully  developed  in 
Phase  II  work. 


ESSENTIAL  TASK  ELEMENTS  -  PHASE  II 

-  Formulate  the  generalized  equations  of  motion  for  the  vessel/fender 
dynamic  interaction  problem  based  on  the  approach  identified  in 
Phase  I  work.  This  approach  will  consider  fender  performance 
algorithms,  local  vessel  stiffness,  dock  mass  and  stiffness,  vessel 
and  berthing  characteristics. 

-  Characterize  vessel  local  hull  or  appendage  stiffness. 

-  Characterize  dock  stiffness  and  mass  characteristics. 

-  Characterize  hydrodynamic  mass  and  damping  for  vessels  considered. 

-  Computer  cod?,  methodology. 

-  Validate  results  against  existing  experimental  data. 

-  Validate  results  against  proposed  test  program. 


II.  DATA  BASE  ACQUISITION  RELATIVE  TO  FENDER-RELATED  INFORMATION 


A  data  base  of  fender  performance  data  and  related  Information  has  been 


established  during  the  initial  efforts  of  this  project.  Some  110  reports, 
papers  and  manufacturers’  catalogs  have  been  accumulated  and  are  included 
as  a  bibliography  with  brief  abstracts  in  Appendix  A. 

The  sources  of  information  relative  to  fender  performance  data  have  been 
designated  by  single  asterisks  (*).  Similarly,  those  sources  relative  to  the 
vessel/fender  response  problem  have  been  designated  with  a  double  asterisk  (**). 

A  search  of  the  available  literature  listed  in  Appendix  A  resulted  in  the 
following  observations  relative  to  fender  performance  data. 

•  The  most  significant  source  of  energy  and  load  deflection  data  is  con¬ 
tained  within  the  catalogs  of  individual  marine  fender  manufacturers. 

•  Fender  performance  data  listed  within  individual  papers  or  reports 
generally  reflect  data  extracted  from  manufacturers’  catalogs. 

•  The  most  common  generic  type  of  fender  for  which  performance  data 
is  available  is  rubber. 

•  A  minimum,  and  in  some  cases  negligible,  amount  of  performance  data 
is  available  for  wood,  gravity,  torsion,  hydropneumatic,  hydraulic 
and  spring  fendering  systems. 

•  The  second  and  third  most  common  fendering .systems  for  which  performance 
data  is  available  is  pneumatic  bag  and  foam-filled  fenders.  There 
appears  to  be  two  primary  manufacturers  of  pneumatic  fenders,  Sumitomo 
and  Yakoharaa,  and  two  primary  manufacturers  of  foam-filled  fenders, 

Seward  and  Samson. 

•  Much  of  the  performance  data  for  large  size  fenders,  especially  pneumatic 
bag  types,  is  extrapolated  and  not  the  result  of  full-scale  testing  due 
to  the  magnitude  of  the  loads  required  in  compression. 


*% 


III.  GENERIC  TYPE  FENDERS  INVESTIGATED 


The  result  of  extensive  literature  search  relative  to  ships  fendering 
systems  concluded  that  fendering  systems  generally  fall  into  two  categories: 
one  considered  commonly  available,  the  other  highly  specialized.  Common 
type  fenders  are  readily  obtainable  from  many  marine  fender  manufacturers 
who  have  performed  extensive  full-scale  fender  tests  relative  to  energy 
absorption  and  load  deflection  data.  These  systems  are  widely  used  for  com¬ 
mercial  and  naval  applications.  Other  systems  considered  highly  specialized 
have  been  determined  to  have  little  performance  data  available  and  have  limited 
or  questionable  practical  field  application.  By  far  the  most  common  fender 
system,  for  which  extensive  performance  data  was  available,  was  rubber  fenders 
followed  by  pneumatic  and  foam-filled  fender  types.  Specialized  torsional, 
hydropneumatic,  gravity  and  hydraulic  fenders  had  little  available  performance 
data  and  relatively  limited  practical  field  application.  The  small  quantity  of 
data  results  available  for  specialized  fender  types  was  concluded  so  specific 
and  unique  to  the  system  investigated  that  generalized  fender  performance  relation 
ships  could  not  be  readily  determined.  Generalized  relationships  require  per¬ 
formance  data  for  variations  of  fender  system  parameters.  Thus,  the  performance 
data  derived  from  a  particular  specialized  fender  test  could  not  be  generalized 
to  describe  the  generic  family-type  action. 

For  the  commonly  available  systems  consisting  of  rubber,  pneumatic  and 
foam-filled  systems,  sufficient  data  was  available  for  variations  of  system 
parameters,  including  basic  geometric  dimensions,  materials,  pressures,  etc. 
in  addition  to  full-scale  test  and  extrapolated  load  and  energy  deflection 
performance  data.  Because  of  the  availability  of  required  information,  rubber 
pneumatic  bag  and  foam-filled  fenders  were  selected  for  detailed  investigation. 
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were  the  following  specific  type  fender  systems 


Included  in  this  investigation 

•  Rubber 
“  Hollow  cylinder  -  transverse  loading 

-  Hollow  cylinder  -  axial  loading 

-  Trapezoidal  -  transverse  loading 

-  Solid  cylinder  -  shear  loading 

-  Hollow  cubic  -  transverse  loading 

-  Hollow  cubic  -  shear  loading 

~  Rotary  donut  -  transverse  loading 

•  Pneumatic 

-  Floating  bag  -  transverse  loading 

-  Air  block  fenders  -  transverse  loading 

-  Air  block  cushions  -  transverse  loading 

•  Foam  filled 

~  Floating  bag  -  transverse  loading 

Table  1  indicates  the  various  manufacturers  of  the  fender  systems  investigated 
in  this  task. 


MANUFACTURER /FENDER  AVAILABILITY 
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IV.  DERIVED  FENDER  PERFORMANCE  RELATIONSHIPS 


APPROACH 

Figure  1  illustrates  typical  energy  and  load/deflection  performance  data 
for  a  rubber  fender  system.  In  this  case  the  data  reflects  the  performance  of 
a  set  of  hollow  cylindrical  fenders  under  transverse  loading.  This  data  is 
typical  of  data  obtained  from  various  manufacturers  and  test  reports.  As  can 
be  seen,  the  energy  and  load  relationships  vary  as  functions  of  geometric  and 
material  properties  between  various  manufacturers. 

In  order  to  characterize  a  particular  type  fender  system,  it  is  necessary 
to  determine  a  relationship  between  the  fender  type  variables  capable  of  con¬ 
densing  or  collapsing  the  particular  performance  relationships  shown  in  Figure  1 
into  more  generalized  ones  which  represent  a  family  of  curves.  An  equation  for 
this  generalized  relationship  can  then  be  determined  as  a  function  of  the  vari¬ 
ables  established.  The  generalized  relationships  are  based  on  the  performance 
data  acquired  from  numerous  sources  and  are  the  basis  for  characterizing  the 
generalized  performance  of  the  fender  type. 

The  resulting  generalized  relationships  derived  will  have  an  associated 
degree  of  dispersion  in  the  relationship  which  reflects  the  choice  of  variables 
selected  to  collapsed  performance  curves  and  the  accuracy  of  the  test  or  extrapo¬ 
lated  performance  data  in  addition  to  the  effects  from  manufacturers’  material 
differences.  The  more  accurate  the  choice  of  collapsing  and  nondiroensionalizing 
variables  and  the  more  accurate  the  available  test  or  extrapolated  data,  the  less 
dispersion  will  be  evident  in  the  relationships  determined. 

The  collapsing  variable  identification  process  is  illustrated  typically 
in  Figure  2  which  correlates  the  energy-absorbing  capacity  of  rubber  cylindrical 
fenders  at  the  rated  50  percent  deflection  to  the  volume  of  material  tested. 


Figure  1:  TYPICAL  PERFORMANCE  DATA 

Load/Deflection  /Energy  Absorption 


Cylindrical  Hanging  Fenders  5  to  24  O.D. 

Load  applied  perpendicular  to  bore.  Test  length:  1  ft 
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Intuitively,  it  is  reasoned  that  a  quantifiable  relationship  between  energy 
capacity  and  volume  exists  but  is  not  evident,  since  energy  storage  per 
unit  volume  may  be  significantly  affected  by  manufacturers’  brand  materials 
and  by  the  variable  ratio  of  cylinder  inner  and  outer  diameters  and  any  other 
pertinent  factors.  However,  the  high  degree  of  linearity  illustrated  by  the 
resulting  relationship  from  four  manufacturers'  data  sources  based  on  some  70 
cylinder  sizes  indicates  that  this  relationship  is  well  defined,  highly  linear 
and  has  very  little  associated  dispersion,  regardless  of  the  material  differences 
and  variation  of  diameter  ratio  generally  existing.  This  identifies  energy 
storage  per  unit  volume  as  a  strong  collapsing  relationship  for  this  fender 
type.  (This  relationship  has  prevailed  for  most  types  of  rubber  fenders 
investigated.)  The  above  correlation  and  significant  condensation  of  particular 
energy  performance  curves  which  occurs  based  on  this  relationship  substantiates 
the  choice  of  the  volumetric  relationship  as  a  significant  relationship  for 
rubber  fender  types.  A  relative  independence  of  the  effects  of  geometry  and 
material  differences  is  implied  by  this  relationship  for  fenders  presently 
available  from  fender  manufacturers.  Figures  3  and  4  illustrate  the  typical 
effect  of  a  valid  and  invalid  choice  of  parameter  selection.  In  Figure  3  it 
is  evident  that  considerable  dispersion  results  from  a  noneffective  choice 
of  condensing  parameters  in  contrast  to  Figure  4  where  the  choice  has  resulted 
in  a  more  well-defined  relationship. 

Once  a  satisfactory  set  of  condensing  or  collapsing  variables  has  been 
formulated  and  determined  to  result  in  a  minimum  of  dispersion  for  the  generic 
relationship  when  plotted,  a  polynomial  regression  analysis  is  performed  to 
determine  the  polynomial  equation  of  order  (n)  which  best  describes  the  energy/ 
deflection  or  load/deflection  generalized  relationship  for  the  fender  type. 
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Figure  3:  NON-EFFECTIVE  PERFORMANCE  CONDENSATION 
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The  equations  which  result  are  considered  characteristic  of  the  performance 
relationships  which  exist  for  the  fender  type  based  on  available  manufacturers' 
data. 

POLYNOMIAL  PERFORMANCE  EQUATIONS  -  GENERAL 

The  polynomial  equations  derived  from  regression  analysis  of  manufacturers' 
data  for  energy  and  load /deflection  are  of  the  following  form: 


rfl  i 

E  *  a  •< 

|E  xix  | 

P  -  T  J  2j  V1  ■ 

l  i-0  J 

L  i-0 

(1) 


where : 


E  *  Fender  energy  absorption 

a  *  Characteristic  fender  volume 

P  *  Fender  reaction  load 

Y  =  Characteristic  fender  area 

X,Y  =  Nondiraensionalized  deflection —  ^ 

(L)  *  characteristic  dimension  of  fender  type 

n  *=  Order  of  polynomial  equation  used  in  regression  analysis 

^i’^i  *  Polynomial  coefficients 


Tables  2  and  3  located  at  the  end  of  this  section,  summarize  the  polynomial 
coefficients  determined  for  energy  absorption  and  load  deflection  described  in 
detail  in  the  following  pages. 
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CYLINDRICAL  FENDERS  -  TRANSVERSELY  LOADED 


Antwerp  Dockt,  Belgium 


Cylindrical  lenders  0  3P0«  IWxTOOP  mm  long 


CYLINDRICAL  FENDERS  -  TRANSVERSE  LOADING 


This  type  fender  is  by  far  the  most  commonly  available  from  various  fender 
manufacturers  as  indicated  in  Table  1.  The  variables  determined  most  effec¬ 
tive  in  condensing  the  test  energy  absorption  performance  data  were  cylinder 
volume  for  energy  absorption,  and  cylinder  inside  diameter  for  nondimensionalizing 
deflection.  Figure  6  indicates  the  generic  energy  absorption  relationship 
which  resulted  from  the  data  sources  considered.  In  this  case,  the  length  of 
cylinder  considered  is  on  a  per  foot  basis.  The  relationship  between  variables 
which  best  fits  the  trend  indicated  in  Figure  7  for  energy  absorption  is: 

E  *  6L  {0.09X  -  5.07X2  +  9.14X3}  103  (2) 

where : 

E  *  Energy  absorption  (ft-lb) 

Dq  *  Outside  diameter  of  the  cylindrical  fender  (ft) 

D^  »  Inside  diameter  of  the  cylindrical  fender  (ft) 

L  *  Length  of  fender  (ft)  (plotted  in  Figure  6  per  foot  length) 

X  *  A/D^  nondimensional 
A  *  Deflection  under  load  (ft) 

6  -  f  O.o2  -  Dl2)  (ft2) 

The  above  equation  is  representative  in  the  range  of  X  <_  1.5. 
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The  load/deflection  relationship  indicated  in  Figure  7  can  be  characterized 
by  the  following  equation. 

P  -  DqL  {105. 76X  -  254. 88X2  +  163. 95X3}  103 
where : 

P  *  Reaction  load  (lb) 

Dq  *=  Cylinder  outside  diameter  (ft) 

X  «=  A/D^  nondimensional 

Di  *  Cylinder  inside  diameter  (ft) 

A  *  Deflection  under  load  (ft) 

L  =  Length  of  fender  (ft)  (plotted  in  Figure  7  per  foot  length) 

The  above  equation  is  representative  in  the  range  of  X  <  1.5. 


Figure  6:  GENERIC  PERFORMANCE  PLOT 

TYPE  RUBBER  CYLINDRICAL 

LOADING:  TRANSVERSE 

PLOT:  ENERGY  ABSORPTION 
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Figure  7:  GENERIC  PERFORMANCE  PLOT 

TYPE  RUBBER  CYLINDRICAL 

LOADING:  TRANSVERSE 

PLOt:  LOAD  /  DEFLECTION 

Cper  foot  length  of  fender) 
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AXIAL  LOADED  CYLINDRICAL  FENDERS 


AXIAL  LOADED  CYLINDRICAL  FENDERS 

Next  to  transversely  loaded  cylinders,  the  axial  loaded  fender  was  the 
most  common  type  for  which  performance  data  was  available.  These  fender 
systems  are  annular  columns  which  compress  and  deflect  as  buckling  columns 
with  added  energy  capability  resulting  from  its  "hoop"  effect.  Being 
circular,  they  have  equal  shear  resistance  for  all  directions  of  transverse 
loading. 

Figure  9  indicates  the  generic  relationship  for  energy  absorption 
determined  for  this  type  fender  system.  The  energy-volume  relationship 
determined  for  transversely  loaded  cylinders  was  equally  effective  for 
axial  fender  types.  In  this  case,  the  nondimensionalizing  characteristic 
dimension  for  deflection  was  the  length  of  the  cylinder.  The  resulting 
algorithm  for  energy  absorption  derived  from  Figure  9  was  determined  to.be: 

E  «  8H  {5.95X  +  51.13X2  +  20.79X3}  103  (4) 

The  corresponding  load/deflection  algorithm  illustrated  in  Figure  10  is: 

P  -  8H  {140.69X+  6.4X2  -  15.65X3}  103  (5) 

where : 

P,  E,  Do,  D^,  B,  A  are  previously  defined. 

X  -  A/H 

H  *  Height  of  fender  (ft)- 


Hie  applicability  of  the  above  equations  is  X  <  0.6. 


Figure  9:  GENERIC  PERFORMANCE  PLOT 


Figure  10:  GENERIC  PERFORMANCE  PLOT 


TYPE 

LOADING 

PLOT 


RUBBER  CYLINDRICAL  FENDERS 
AXIAL  COMPRESSION 
LOAD  /  DEFLECTION 


d  UNIROYAL 
A  U.S.  RUBBER 
o  GOODYEAR 


\ 

i. 


HOLLOW  CUBIC  RUBBER  FENDERS  -  SHEAR  LOADED 
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■BGLTSte)  1/2  UNC  -2  A  x  1-3/4  LG 
CORROSION  PROTECTION  REQD. 
SAE  GRADE  5  ' 

■STANDARD  WASHER 


.  SHEAR  FENDER 


DOCK  STRUCTURE  TYR  (FIXED) 


HOLLOW  CUBIC  ROBBER  FENDERS  -  SHEAR  LOADED 

Shear  fenders  have  the  ability  to  stretch  in  four  shear  directions  in 
addition  to  withstanding  large  compression  loads  and  limited  tension  loads. 
This  feature  allows  for  wale  movement  away  from,  into  or  tangential  to  docks 
as  vessels  berth.  In  compression,  the  fender  can  be  used  alone  or  in  tandem, 
bolted  between  a  wall  and  whale.  Tension  and  compression  loading  allow  the 
shear  fender  to  support  the  wale.  This  is  illustrated  in  the  previous  figure 
(Figure  11).  Although  this  type  fender  system  is  simple  and  effective,  it  is 
not  commonly  available  from  fender  manufacturers  as  indicated  in  Table  1. 

The  data  used  to  determine  the  appropriate  performance  algorithms  has 
been  selected  from  a  single  source  but  reflects  the  eight  different  size 
fenders  available.  These  fenders  are  characterized  by  a  cylindrical  bore 
running  lengthwise  in  the  direction  of  shear  loading. 

Figure  12  indicates  the  characteristic  energy  absorption  curve  deter¬ 
mined  by  correlating  the  fender  energy  volume  and  the  normalized  deflection 
relationship.  In  this  case  the  characteristic  length  was  the  height  of  the 
shear  fender  normal  to  the  direction  of  shear.  The  representative  equation 
for  energy  absorption  was  determined  to  be: 

E  -  6H  {2.63X  -  5.39X2  +  10.62X3  -  3.44X4}  103  (6) 

The  corresponding  load/deflection  relationship  illustrated  in  Figure  13  is: 

F  *  8  {21. 6X  -  21. 9X2  +  19.76X3  -  5.06X4}  103  (7) 

where : 

P  *  Reaction  load  (lb) 

B  -  *bMBd 2  <ft2) 
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W.  *  Width  of  fender  base  (ft) 

D 

B^  *  Bore  diameter  (ft) 

H  *  Height  of  shear  fender  (ft) 

X  *  A/H  (nondimensional) 

E  *  Energy  absorption  (ft-lb) 

The  above  equations  are  representative  for  X  £  1.9 


Figure  12:  GENERIC  PERFORMANCE  PLOT 


TYPE:  HOLLOW  RUBBER  CUBIC 
-LOADING:  SHEAR 


PLOT:  ENERGY  /  DEFLECTION 
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Figure  13:  GENERIC  PERFORMANCE  PLOT 
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HOLLOW  CUBIC  RUBBER  FENDERS  -  TRANSVERSELY  LOADED 


Figure  14:  HOLLOW  CUBIC  COMPRESSION  FENDER  INSTALLATION 
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HOLLOW  CUBIC  RUBBER  FENDERS  -  TRANSVERSELY  LOADED 

This  fendering  system  generally  offers  larger  energy  absorption  and 
larger  reaction  loads  compared  to  similar  sized  cylindrical  fenders  although 
it  does  not  exhibit  the  typical  buckling  phenomenon.  Their  use  is  generally 
between  wood  walers  and  concrete  piers  or  draped  as  indicated  in  Figure  14. 

The  parameters  determined  significant  in  collapsing  the  energy  and 
reaction  load  relationships  were  determined  to  be  fender  volume  and  fender 
height  in  the  loaded  direction.  Fender  deflections  were  normalized  by  the 
characteristic  fender  height. 

Figure  15  indicates  the  resulting  generic  relationship  for  energy 
absorption.  This  curve  can  be  defined  by  the  following  equation. 

E  =  HWbL  {21. IX  -  74. IX2  +  208. 8X3}  103  .  -  (8) 

The  corresponding  load/deflection  relationship  illustrated  in  Figure  16  is: 

P  -  HL  {178. 7X  -  702. 8X2  +  1600. 9X3}  103  (9) 

where : 

E  *  Energy  absorption  (ft-lb) 

P  *  Reaction  load  (lb) 

H  *=  Height  of  cubic  in  direction  of  loading  (ft) 

W,  =  Base  width  of  cubic  (ft) 

D 

L  *  Length  of  fender  (ft) 

X  *  A/H  (nondimens ional ) ,  X  £  0.65 


E  /WUH  x  10-3  CFT-LB/FT3) 


Figure  15:  GENERIC  PERFORMANCE  PLOT 


TYPE:  HOLLOW  RUBBER  CUBIC 

LOADING:  TRANSVERSE 

PLOT:  ENERGY  ABSORPTION 

(per  foot  length  of  fender) 
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Figure  16:  GENERIC  PERFORMANCE  PLOT 


TYPE:  HOLLOW  RUBBER  CUBIC 

LOADING:  TRANSVERSE 

PLOT  :  LOAD/  DEFLECTION 

(PER  FOOT  LENGTH  OF  FENDER) 
O  GOODYEAR 
□  UNIROYAL 
A  JOHNSON  RUBBER 
o  U.  S.  RUBBER 


TRAPEZOIDAL  RUBBER  FENDERS  -  TRANSVERSELY  LOADED 


TRAPEZOIDAL  RUBBER  FENDERS  -  TRANSVERSELY  LOADED 

Trapezoidal  rubber  fenders  employ  two  mechanisms  in  the  absorption  of 
energy,  these  are:  direct  compression  and  buckling.  They  are  generally 
mounted  directly  to  open-faced  structures, or  they  can  be  used  in  combination 
with  timbering. 

Figure  28  indicates  the  generic  load/deflection  curve  which  is 
characterized  by  the  region  of  buckling  generally  occurring  at  approximately 
30  percent  deflection. 

Figure  19  illustrates  the  generic  energy/def lection  curve  for  this  type 
fender.  For  this  case  the  characteristic  height  of  the  fender  in  the  direction 
of  loading  was  determined  to  be  the  significant  condensing  parameter  for  energy 
absorption  and  load/deflection. 

The  following  equations  are  representative  of  tradezoiaal  fender  _ 
performance: 

E  «=  HLW  {0.57X  +  36.55X2  -  56.55X3  +  40.37X4}  103  (10) 

D 

The  corresponding  load/deflection  relationship  is: 

P  =  HL  (105. 82X  -  207. 06X2  -  48.24X3  +  423. 72X4}  103  (11) 

where: 

P  *  Reaction  load  (lb) 

*  Fender  base  width  (ft) 

H  *  Fender  height  in  direction  of  load  (ft) 

L  *  Fender  length  (ft) 

E  *  Energy  absorption  (lb-ft) 

X  *  A/H  (nondimensional,  X  <  0.53) 


Figure  18:  GENERIC  PERFORMANCE  PLOT  © 


TYPE:  TRAPEZOIDAL  RUBBER 
CLOSED  BASE 

LOADING:  TRANSVERSE  COMPRESSION 

PLOT:  LOAD  /  DEFLECTION 
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Figure  19:  GENERIC  PERFORMANCE  PLOT 


TYPE:  TRAPEZOIDAL  RUBBER  -  CLOSED  BASE 
LOADING:  TRANSVERSE  COMPRESSION 


PLOT:  ENERGY/ DEFLECTION 

(per  foot  length  of  fender) 


A/H 

(NON  -  DIMENSIONAL) 
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SOLID  CYLINDRICAL  SHEAR  FENDERS 

In  shear  fender  installations,  as  in  Figure  20,  wooden  fendering  is  fitted 
over  by  means  of  supports  on  metal  plates.  The  fenders  are  then  mounted  on 

i* 

brackets  secured  to  the  quay.  When  the  wood  fendering  is  compressed,  the  shear 
fenders  are  loaded  into  shear.  Since  the  shear  modulus  of  rubber  is  only  a 
third  of  its  modulus  of  elasticity,  reaction  forces  are  kept  low  for  this  type 
configuration. 

Figure  21  indicates  the  generic  energy  absorption  relationship  which 
characterizes  this  type  fender.  In  the  figure  fender  energy  absorption 
curves  have  been  condensed  by  the  volume  of  the  fender.  The  deflection 
under  load  has  been  nondimensionalized  by  the  fender  height  normal  to  the 
loading. 

The  energy  equation  which  characterizes  this  relationship  is: 


E  *  6H  {0.54X  +  8.79X2}  103 


(12) 


The  corresponding  load/deflection  relationship  illustrated  in  Figure  22  is: 


P  =  B  {22.77X  +  1.14X2  -  1.43X3}  103 


where : 
P 

B 


Reaction  load  (lb 
2 


rD 


(ft2) 


\ 

J 


X  *=  A/H,  X  <  1.0 

Dq  *  Diameter  of  cylinder  (ft) 

H  *  Height  of  fender  (ft) 

A  *  Fender  deflection  under  load 
E  *  Energy  absorption  (ft-lb) 


(13) 


44 


0(0 


Figure  21:  .GENERIC  PERFORMANCE  PLOT 


Figure  22:  GENERIC  PERFORMANCE  PLOT 


TYPE:  SOLID  RUBBER  CYLINDER 
LOADING:  SHEAR 
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o  VREDESTEIN 
□  REGAL  ■ 


ROTARY  DONUT  FENDERS  -  TRANSVERSE  COMPRESSION 


ROTARY  DONUT  FENDERS  -  TRANSVERSE  COMPRESSION 

Rotary  fenders  consist  of  hollow  section  rubber  wheels  which  are  mounted 
on  a  control  axis  that  allows  them  to  rotate  freely  when  ships  horizontal  shear 
forces  are  applied.  This  type  fender  is  available  in  multiple-wheeled  configu¬ 
rations  in  a  variety  of  wheel  size  diameters.  The  hollow  wheeled  section 
essentially  absorbs  its  energy  in  material  compression  and  exhibits  signifi¬ 
cant  absorption  compared  to  the  reaction  loads  developed. 

Figures  24  and  25  illustrate  the  energy  absorption  and  reaction  load 
deflection  curves  derived  for  single-,  double-  and  triple-wheeled  fender 
configurations.  In  these  figures  the  energy  absorption  relationships  have 
been  modified  by  characteristic  dimensions  of  fender  inner  and  outer  diameters, 
the  number  of  donuts  per  axial  and  the  width  of  the  donut  base.  For  the 
reaction  load  relationship:  the  number  of  donuts,  outer  diameter  and  base 
width  were  significant.  For  both  relationships  deflection  was  normalized  by 
the  characteristic  depth  of  the  donut  tire. 

The  generic  relationships  derived  for  energy  absorption  and  load/deflection 
and  indicated  in  Figures  24  and  25  are: 

E  -  NB  {5.51X  -  21.31X2  +  28.05X3}  103  (14) 

The  corresponding  load  deflection  curve  determined  was: 

P  =  ND  W,  {-0.45X  +  67.32X2  -  189. 6X3  +  188.46XA)  103  (15) 

o  b 

where: 


N  =  Number  of  donuts  per  axial 


Figure  24:  GENERIC  PERFORMANCE  PLOT 


TYPE:  ROTARY  DONUT 
LOADING:  TRANSVERSE  COMPRESSION 
PLOT:  ENERGY  ABSORPTION 
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Figure  25:  GENERIC  PERFORMANCE  PLOT 

TYPE  ROTARY  DONUT 
LOADING:  TRANSVERSE  COMPRESSION  ■ 

PLOT:  LOAD  /  DEFLECTION 
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PNEUMATIC  RUBBER  FENDERS  -  FLOATING  TYPE 


PNEUMATIC  RUBBER  FENDERS  -  FLOATING  TYPE 

Floating  pneumatic  fenders  utilize  the  compressive  elasticity  of  air 
to  support  loads.  For  this  reason  performance  deterioration  due  to  fatigue 
is  absent.  For  realistic  oblique  ship  loading,  pneumatic  type  fenders  do  not 
suffer  significant  loss  of  energy  absorbing  capacity  as  do  solid  rubber  fenders. 
For  rough  weather  mooring,  this  type  fendering  system  exhibits  much  less  damage 
due  to  the  fact  that  maximum  reaction  forces  under  combined  shear  and  compres¬ 
sion  increase  slowly  and  sustain  large  allowable  deflections.  Under  excessive 
loads  these  fenders  do  not  result  in  excessive  reaction  loads  as  do  solid  or 
bottomed  out  rubber  fenders. 

Figures  27  and  28  indicate  the  results  for  32  different  size  pneumatic 
fenders  investigated.  These  fenders  ranged  in  pressures  from  4.3  to  11.4  psi 
internally.  Figure  27  illustrates  the  characteristic  energy  absorbing  relation¬ 
ship  resulting  from  condensing  the  plot  of  energy  absorption  by  the  relationship 
pl/1.4  ld^  the  pressure,  length  and  diameter  characteristic  of  the  fender. 

This  quantity  is  plotted  against  the  deflection  normalized  by  the  diameter  of 
the  cylinder  bag. 

The  following  relationship  was  determined  representative  of  energy  absorp¬ 
tion  for  pneumatic  fenders. 

E  **  BD  {0.82X  -  2.54X2  +  17.94X3}  103  (16) 

o 

The  corresponding  load/deflection  relationship  illustrated  in  Figure  28  is: 

P  -  6  {5.19X  +  39.95X2  -  77.02X3  +  149. 09X4}  103  (17) 


where : 

6 


pl/1.4 


LD 

o 
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Internal  pressure  (psf) 
Length  of  fender  (ft) 
Fender  diameter  (ft) 

A/D  ,  X  <  0.55 

O  ““ 

Fender  deflection  (ft) 
Energy  absorption  (lb-ft) 


Figure  27:  GENERIC  PERFORMANCE  PLOT 


TYPE:  RUBBER  PNEUMATIC  *  FLOATING  BAG 
LOADING:  TRANSVERSE  COMPRESSION 
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Figure  28:  GENERIC  PERFORMANCE  PLOT 

TYPE:  RUBBER  PNEUMATIC  -  FLOATING  TYPE 
LOADING:  TRANSVERSE  COMPRESSION  £ 

PLOT:  LOAD  /  DEFLECTION 
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RUBBER  PNEUMATIC  -  AIR  BLOCK  FENDER 


RUBBER  PNEUMATIC  -  AIR  BLOCK  FENDER 

Air  block  fenders  are  pneumatic,  axially  loaded  fenders  which  can  be 
bolted  to  docks  and  applicable  when  floating  bag  types  cannot  be  used. 

They  offer  all  the  performance  advantages  that  pneumatic  bag  types  generally 
exhibit. 

The  characteristic  performance  curves  illustrated  in  Figures  30  and  31 
were  determined  by  nondimensionalizing  the  energy  absorption  and  load/deflection 
curves  by  P,  H,  and  D,  the  characteristic  pressure,  height  and  diameter  of  the 
block  fender.  The  energy  absorption  and  load/defelction  relationships  were  then 
plotted  against  nondimensional  deflection  A/H,  the  percentage  fender  height. 

The  resulting  relationships  are  based  on  the  investigations  of  13  fender  sizes 
at  14.2  psl.  Since  this  type  fender  was  available  in  only  one  pressure  size, 
the  pressure  variable  was  considered  similar  to  the  relationship  determined 
for  floating  bag  types. 

The  following  relationship  is  representative  of  the  energy  absorption  of 
air  block  fenders  illustrated  in  Figure  30: 

E  =  8D  {2.58X  +  9.73X2  -  13.40X3  +  40.09X4}  103  (18) 

o 

The  corresponding  load/deflection  relationship  illustrated  in  Figure  31  is 
represented  by: 

P  -  6  (43.96X  -  8.77X2  -  62.48X3  +  256. 23X4}  103  (19) 

where: 

E  «  Energy  absorption 

P  *  Reaction  load 
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6  -  P1/1-*BD0 

H  *  Fender  height 

D  *  Fender  diameter 
o 

p  *  Internal  pressure 
X  =  A/H  (nondimens ional)  X  £  0.6 
A  *  Fender  deflection 

The  above  relationships  are  valid  for  any  set  of  consistent  units. 
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Figure  30:  GENERIC  PERFORMANCE  PLOT 


TYPE:  RUBBER  PNEUMATIC  -  AIR  BLOCK  .FENDERS 
LOADING:  AXIAL  COMPRESSION 
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Figure  31:  GENERIC  PERFORMANCE  PLOT 

TYPE:  RUBBER  PNEUMATIC  -  AIR  BLOCK  FENDERS 
LOADING:  AXIAL  COMPRESSION 
PLOT:  LOAD  /  DEFLECTION 


O  YOKOHAMA 


!  RUBBER  PNEUMATIC  -  AIR  BLOCK  CUSHIONS 

* 


i 

l 

,  t 

$ 


65 


RUBBER  PNEUMATIC  -  AIR  BLOCK  CUSHIONS 


This  system  is  similar  to  the  air  block  fenders  except  that  it  is 
rectangular  in  shape.  It  is  mounted  on  a  steel  backing  plate  which  can 
be  bolted  to  docks  or  semisubmersible  drill  rig  legs.  Although  this  fendering 
system  appeared  to  have  a  significant  number  of  merits  typically  associated 
with  pneumatic  systems,  they  were  only  available  from  one  fendering  manu¬ 
facturer. 

The  data  represented  in  Figures  33  and  34  are  for  only  two  fender 
lengths  at  the  same  internal  pressure.  In  these  figures  the  energy  and 
load  curves  have  been  normalized  by  the  characteristic  pressure,  length 
and  base  width  dimensions,  while  the  deflection  has  been  normalized  by 
the  cushion  height  in  the  direction  of  loading. 

The  derived  relationship  which  best  fits  the  condensed  data  for  energy 
was  determined  to  be: 

E  «  BH  {-0.12X  +  7.46X2  -  12.71X3  +  14.77X4}  103  (20) 

The  corresponding  load  deflection  relationship  is: 

P  »  B  (9.22X  -  4.16X2  +  5.10X3  +  29.55X4}  103  (21) 

where : 

E  =  Energy  absorption 

P  *  Reaction  load 
B  «  p1/3-4  «  L 

D 

p  *  Internal  pressure 

W,  *  Base  width  of  cushion 
b 
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L  *  Cushion  length 

H  *  Cushion  height  in  direction  of  load 
X  >=  A/H,  X  <  0.6 
A  *  Fender  deflection 


The  above  relationships  are  valid  for  any  set  of  consistent  units. 


Figure  33:  GENERIC  PERFORMANCE  PLOT 
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Figure  34:  GENERIC  PERFORMANCE  PLOT 
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FOAM-FILLED  RUBBER  FENDERS 

This  system  resembles  the  pneumatic  bag  type  fenders  except  internally 
they  are  completely  filled  with  resilient,  closed-cell  foam.  They  typically 
have  greater  energy  absorption  and  less  reaction  force  than  pneumatic  fenders 
of  equal  size.  They  are  generally  lighter  than  pneumatic  fenders  of  equal 
capacity  and  cannot  explode  or  sink  if  punctured. 

Figures  36  and  37  indicate  the  results  of  investigating  two  primary 
sources  of  performance  data  for  this  type  fender  system.  For  these  curves 
the  parameters  found  to  condense  the  performance  relationship  were  fender 
length  and  diameter. 

The  energy  absorption  relationship  found  characteristic  of  Figure  36  was 
determined  to  be: 

E  «  BD  {0.27X  -  1.03X2  +  6.43X3  -  A.69X4}  103  '  (22) 

o 

The  corresponding  load  deflection  relationship  illustrated  in  Figure  37  is: 

P  =  6  {1.77X  +  6.25X2  -  13.81X3  +  16.32X4}  103  (23) 


where : 

E  = 
P  * 

B  * 

D  - 
o 

L  = 
X  * 
A  * 


Energy  absorption  (ft-lb) 

Reaction  load  (lb) 

D  L 
o 

Fender  diameter  (ft) 

Fender  length  (ft) 

A/D  (ncndimensional) 
o 

Fender  deflection  (ft) 
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Figure  36:  GENERIC  PERFORMANCE  PLOT 
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SUMMARY  TABLES  FOR  GENERIC  FENDER  ALGORITHMS 
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V.  RANKING  OF  FENDER  SYSTEM  MECHANISMS 


Commonly  available  fenders  from  manufacturers  of  fender  systems  operate 
on  the  basis  of  one  or  more  mechanisms  which  determine  the  way  in  which  the 
fender  stores  energy  and  deflects  under  loading.  These  mechanisms  are  commonly: 

a.  Axial  compression 

b.  Transverse  compression 

c.  Transverse  shear 

d.  Pneumatic  bag  compression 

e.  Foam-filled  bag  compression 

These  different  mechanisms  result  in  considerable  differences  in  the  basic 
performance  characteristics  of  the  individual  fender  types.  The  ranking  of 
energy-absorbing  mechanisms  takes  on  a  significantly  different  importance 
depending  on  the  measure  of  merit  or  goal  which  is  established  for  the  ranking 
process.  Since  a  designer  is  concerned  with  many  variables  such  as  fender 
energy  absorption,  reaction  load,  deflection,  relative  system  costs,  system 
durability,  etc. ,  the  ranking  of  fender  mechanisms  will  vary  in  accordance 
with  his  selected  criteria.  For  purposes  of  this  discussion,  only  two  measures 
of  merit  are  considered:  energy  absorption  capability  and  reaction  load  as  a 
function  of  deflection.  These  measures  of  merit  are  generally  diametrically 
opposed. From  a  design  point  of  view,  one  would  like  maximum  energy  absorption 
with  minimum  reaction  load  generation  for  a  given  deflection.  In  ranking  the 
candidate  mechanisms,  the  first  approach  considers  which  mechanisms  absorbed 
the  most  energy  for  a  given  deflection  with  reaction  load  not  a  factor.  The 
second  viewpoint  considers  which  mechanisms  resulted  in  the  least  reaction 
load  for  a  given  deflection  not  considering  energy  absorption. 
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VI.  PRELIMINARY  FORMULATION  OF  THE  FENDER/VESSEL 
INTERACTION  RESPONSE  PROBLEM 


The  references  included  in  Appendix  A  relative  to  the  ship/fender 
response  problem  have  been  designated  by  (**) .  These  references  approach 
the  dynamic  response  problem  in  various  ways  and  to  various  depths.  Of  the 
references  cited  in  Appendix  A,  the  "Dynamic  Response  of  the  Ship  and  the 
Berthing  Fender  System  after  Impact,"  (37)  included  as  Appendix  B  for 
ready  reference,  was  considered  the  most  appropriate  for  further  development. 

The  response  problem  formulation  appears  generalized  enough  to  be 
adapted  to  include  the  generic  fender  algorithms  preliminarily  developed 
in  Phase  I  work  and  hull,  dock  and  berthing  characterizations. 


The  essential  task  steps  envisioned  for  Phase  II  efforts  would  include: 

•  Formulate  the  generalized  equations  of  motion  for  the  vessel/fender 
dynamic  interaction  problem  based  on  the  approach  identified  in 
Phase  I  work.  This  approach  will  consider  fender  performance 
algorithms,  local  vessel  stiffness,  dock  mass  and  stiffness,  vessel 
and  berthing  characteristics. 

•  Characterize  vessel  local  hull  or  appendage  stiffness. 

•  Characterize  dock  stiffness  and  mass  characteristics. 

•  Characterize  hydrodynamic  mass  and  damping  for  vessels  considered. 

•  Computer  code  methodology. 

•  Validate  results  against  existing  experimental  data. 

•  Validate  results  against  proposed  test  program. 

It  is  envisioned  that  the  first  two  task  elements  above  would  be  based  on 
Phase  I  results,  references  (37,  (32)  and  (33),  the  basic  methodology  for 
the  dynamic  problem  and  studies  related  to  local  hull  stiffnesses.  Task  3 
will  be  approached  through  a  representative  dock  characterization  for  the 


ship  selected  initially  as  part  of  the  response  problem.  The  hydrodynamic 
mass  and  damping  characteristics  for  this  vessel  would  be  investigated  using 
references  (37),  (20)  and  (75)  in  addition  to  other  relevant  sources  of  mass 
and  damping  information. 

It  is  assumed  that  computer  coding  of  the  dynamic  equations  of  motion 
and  their  solution  will  require  the  significant  Phase  II  effort.  An  initial 
validation  effort  will  include  correlation  between  program  results  and  any 
known  test  results  for  which  comparisons  can  be  made.  These  will  consider 
the  results  in  references  (8),  (62),  (21),  (26)  and  (3)  but  not  be  limited 
to  those  references. 

Actual  test  programs  to  be  developed  as  part  of  task  7  would  consider 

validation  of  fender  algorithms  for  large  size  generic  fenders  via  static 

»  ** 

or  model  testing  since  most  data  issued  by  manufacturers  is  based  on  extrapolation 
of  small  scale  test  data.  This  would  be  further  developed  as  part  of  Phase  II 
efforts.  In  addition,  test  programs  could  include  validation  of  response 
program  results  through  small  scale  model  testing.  This  also  would  be 
developed  further  into  Phase  II  efforts. 
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for  the  solutions  of  the  dynamic  equations  are  included.  These  are  comprised  of 
the  virtual  mass  of  the  ship,  in  both  translational  and  rotational  motion,  in 
addition  to  the  time  interval  required  for  the  solutions  of  the  motion  equations 
by  numerical  integration  methods. 
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by  berthing  and  mooring  ships  have  significant  effects  on  energy-absorption 
capabilities  of  conventional  timber  fender  piles.  This  paper  presents 

•  significant  findings  of  on-site  investigations  and  subsequent  technical  and 
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Second  Series."  CEL,  March  1967.  This  report  summarizes  the  results  of  in- 
service  tests  of  two  experimental  hydraulic-pneumatic  floating  fenders  in  the 
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unsheltered  harbors,  standard  greenheart  timber  pile  with  rubber  bearing  block  at 
deck  level;  and  (3)  for  dock  corners  generally,  the  Raykin  (rubber-in-shaer) 
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extensive  research  and  development  program  executed  at  the  U.S.  Naval  Civil 
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systems;  advantage  and  disadvantages  of  various  fender  systems,  cost- 
effectiveness  and  design  procedures  for  different  marine-  environment  and  exposure 
conditions.  The  energy  absorption  characteristics,  berthing  velocity,  and 

virtual  mass  of  ship  are  discussed  in  detail. 
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i  tanker  to  a  jetty  is  presented.  The  method,  which  results  in  the  perdiction  of 
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International  Corporation  for  "Seibu"  Dock  Fender  in  North  America.  In  this  paper 
he  reviews  three  systems,  the  ship  (docking  operation),  the  fender  (design  and 
spacing),  and  the  berthing  structure. 


|  58.  Powell,  R.G.,  Carle,  R.B.,  “The  Use  of  Hydraulic  Cushioning  in  the  Docking  of 
i  Super  Tankers."  Offshore  Technology  Conference,  May  1972. 

i 

1  1*59.  Ouinn,  A.,  "Design  and  Construction  of  Ports  and  Marine  Structures."  McGraw- 
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ABSTRACT 

ys 

This  paper  describes  a  mothodfof  analysis  for 
evaluating  the  portion  of  ship  kinetic  energy  and 
impact  force  transmitted  to  a  berthing  structure 
provided  with  fenders  which  have  linear  or  non¬ 
linear  spring  constants^  In  the  analysis,  present-  / 
ed  herein, /the  dynamic  responses  of  the  ship,  CL' 
fender  and  Berthing  structure,  after  impact,  are  ^ 
considered,  and  derived  equations  for  the  selec¬ 
tion  of  different  parameters  needed  for  the  solu- 
tions  of  the  dynamic  equation^  are  included. 
These  are  comprised  of  the  virtual  mass  of  the 
ship,  in  both  translational  and  rotational  motion, 
in  addition  to  .the  time  interval  required  for  the 
solutions  of  the  motion  equations  by  numerical 
integration  methods. 


1.  INTRODUCTION 

Since  the  size  of  ships,  particularly  tankers, 
has  increased  in  recent  years,  the  design  of  off¬ 
shore  berthing  structures  has  become  more  im¬ 
portant.  One  of  the  prime  difficulties  facing  de¬ 
signers  is  the  evaluation  of  the  portion  of  .thip 
kinetic  energy  and  the  impact  force  transmitted 
to  each  of  the  fenders  and  the  berthing  struc¬ 
ture,  especially  structures  provided  with  rubber- 
like  fenders. 

( 1 )  The  Kinetic  Energy  of  the  Berthing  Ship 

When  a  ship  is  approaching  the  berth  with 
both  translational  and  rotational  motion,  its  kine¬ 
tic  energy  is  given  by  the  following  equation; 

J . (1) 
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where  Mi—' virtual  .mass  of  the  ship, 

le  ^virtual  moment  of  inertia  about  the 
vertical  axis  through  the  ship's  cen¬ 
ter  of  gravity, 

V»  as  velocity  of  translation, 
tut  =angular  velocity. 

(2)  The  Effective  Energy  for  Fender  System 
Design 

During  berthing  the  kinetic  energy  of  the  ship 
may  be  dissipated  in  several  ways,  among  which 
are  the  following: 

i)  Elastic  deformation  of  the  structure  and 
fender. 

ii)  Swinging  of  the  ship  due  to  yawing  mo¬ 
tion. 

iii)  Heeling  of  the  ship  due  to  rolling  motion. 

iv)  Elastic  deformation  of  the  ship's  hull. 

'v)  Piling  of  the  water  trapped  between  the 
ship's  hull  and  the  face  of  the  berthing 
structure.  (This  occurs  in  the  case  of  a 
long  closed  structure.) 

Designers  who  are  involved  with  marine  struc¬ 
tures  design  are  interested  in  the  portion  of 
energy  indicated  by  i)  which  is  called  the  effec¬ 
tive  energy  {£«).  The  problem  of  determining 
the  effective  energy  has  been  treated  analytically 
by  several  investigators.  Miehalos”*”  treated  the 
problem  as  one  which  had  a  single-degree-of- 
freedom  dynamic  motion,  and  considered  the 
theory  of  elastic,  impact  in  his  analysis.  The 
judgement  of  others”*”,  including  the  authors, 
is  that  the  ship’s  dynamic  impact  on  the  struc¬ 
ture  can  be  considered  as  a  plastic  impact,  where, 
upon  impact,  both  the  ship  and  the  fender  sys¬ 
tem  move  together  as  one  combined  mass.  Vasco 
Costa”  has  derived  a  dynamic  equation  for  esti¬ 
mating  the  effective  energy  in  which  only  the 
portion  of  the  energy  dissipated  by  the  yawing 
motion  of  the  ship  was  considered.  The  rolling 
motion  and  the  influence  of  the  fender  system 
dynamic  response  were  ignored.  Hayashi  and 
Shirai”  have  dealt  with  the  problem  as  one  which 
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has  three-degrees-of-freedora  dynamic  motion; 
swaying,  yawing,  and  rolling  motion.  The  equa- 
tions  are  valid  i)  for  structures  which  are  pro¬ 
vided  with  linear  spring  constants,  such  as  steel 
spring-like  fenders,  ii)  for  only  one  case  of  ap¬ 
proaching  mode  of  berthing,  in  which  the  vector 
of  the  approaching  velocity  is  perpendicular,  to 
the  arm  connecting  the  ship's  center  gravity  with 
the  point  of  contact. 

Besides,  the  dynamic  response  of  the  fender 
system  was  ignored.  The  empirical  equation  for 
determining  the  effective  energy  (£,)  is  also  used 
for  design  purposes  and  is  of  the  form 

E.-CE,  . (2) 

where  E»  represents  the  approaching  ship's  kine¬ 
tic  energy  and  C  is  the  reduction  or  dispersion 
factor.  Pages'0  suggests  the  following  equation 
for  determining  for  C; 

C=l/(l+16o«)  . (3) 

where  e—d/L.  L  represents  the  ship’s  total  length 
and  d  represents  the  distance  between  the  ship's 
center  of  gravity  and  the  point  of  contact,  meas¬ 
ured  parallel  to  the  berthing  face.  Other  design¬ 
ers  have  selected  a  value  of  C  which  varies  from 
0.2  to  1.0  depending  on  several  factors,  such  as 
the  mode  of  berthing  operation,  local  hull  de¬ 
formation,  structure  type,  etc.1'*0-0-10 

From  published  information  it  became  clear 
that  the  portion  of  the  energy  transmitted  sepa¬ 
rately  to  the  berthing  structure  and  fender  is 


still  a  problem,  especially  in  the  case  of  fenders 
with  non-linear  spring  constants,  and  rubber-like 
fenders,  which  are  in  this  class,  are  being  used 
extensively  due  to  their  large  energy-absorbing 
characteristics. 

The  authors  have  presented  a  method  of  analy¬ 
sis  based  on  the  dynamic  behavior  of  the  system, 
after  collision,  to  evaluate  the  impact  load  and 
the  portion  of  the  ship's  kinetic  energy  trans¬ 
mitted  to  the  berthing  structure  and  fenders 
which  includes  fenders  that  have  both  linear  and 
non-linear  spring  constants.  Also,  to  evaluate 
the  portion  of  the  energy  dissipated  in  the  swing¬ 
ing  and  rolling  of  the  ship  after  impact. 

2.  DYNAMIC  RESPONSE  OF  THE  SHIP 
AND  THE  BERTHING  STRUCTURE 
AFTER  THE  FIRST  IMPACT 

( 1 )  General  Mode  of  Berthing 

When  the  ship  is  approaching  the  berth  under 
its  own  power,  it  is  angled  in  to  make  the  first 
contact  with  the  fender  system  at  a  point  near 
its  bow  or  stern.  This  point  of  contact  is  always 
located  in  a  horizontal  plane  higher  than  that 
passing  through  the  ship's  center  of  gravity. 
During  this  mode  of  berthing,  the  ship  will 
undergo  dynamic  motion  which  has  three-degree” 
of -freedom,  namely;  swaying,  yawing,  and  roll¬ 
ing.  The  other  motions,  heaving,  pitching  and 


Fig.  1  Behaviour  of  ship  after  contact. 
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cursing  are  of  little  consequence  in  energy  dis¬ 
sipation  and  may  be  neglected.  • 


( 2 )  Equations  of  Motion 

In  the  following  analysis  it  is  assumed  that  no 
sliding  contact  is  made  along  the  fender's  surface. 

Consider  the  motion  of  the  ship,  as  a  free  dody 
under  the  action  of  the  load  P  acting  at  the  point 
of  contact  C,  we  have; 

At  any  time  t,  after  the  ship  came  in  contact 
with  the  fender  system,  Fig.  1,  its  center  of 
gravity  will  sway  in  the  direction  of  the  acting 
load,  to  the  position  C.  Then  the  ship  will  swing 
about  its  vertical  axis  passing  through  C,  an 
angle  0  due  to  the  yawing  motion  and  finally  it 
will  roll  about  its  horizontal  longitudinal  axis  an 
angle  p.  Denote  the  coordinates  of  the  final 
positions  of  G  and  C,  with  respect  to  the  axes 
X  and  Y,  by  (Xg,  Yg)  and  {Xo,  Yc)  respectively. 
The  X  and  Y  axes  are  taken  parallel  and  nor¬ 
mal  to  the  surface  of  berthing  respectively. 
From  the  figure,  the  relation  connecting  the 
fender  system  motion  at  C  with  the  motion  of 
G  is;  • 

Xg— Xg+(r!+Hp)  cos  (f+a)  | 
?c=¥o-[r!+Hi)s\n{r+a)  j  . (<) 

the  second  order  terms  in  the  above  equation, 
d1,  are  neglected. 

Consider  the  dynamic  equilibrium  of  G,  the 
following  equations  of  motion  will  hold; 


i)  SWAY 

MtXg—  —Px—Rwx  ■. .  • 
M,?g=-Py-Rwy 

ii)  YAWING 
/j_jd=/V-r-sin  (r+o) 

-ParT-cos  (y+o)-M 

iii)  ROLLING 

/i->(i=(PfCOs  ip— P*-sin  <p)H 

-W-Xrt-N, 


(5) 


Px,  Pr  and  Rwx,  Rwy  are  the  components  of  the 
render  reaction  and  the  water  resistance,  after 
:he  ship  came  in  contact  with  the  fender  system, 
•espectively.  Ni  and  Nt  denote  the  water  resis- 
a  nee  to  yawing  ar.d  rolling  motions  respectively. 
However,  as  the  time  of  contact  is  very  small, 
water  resistance,  is  safely  allowed.  Water  rests- 
ince  is  effective  in  the  time  between  the  first 
ind  the  second  impact,  this  will  be  discussed  in 
letails  in  next  paper. 

As  for  the  berthing  structure  response,  the 
ollowing  eq.  will  hold; 


MiPi*=Py-KyY,-flyY,  j  . (  ' 

In  the  previous  equations  the  value  of  p  is  evalu¬ 
ated  from  the  given  load — deflection  relation  (P- 
D)  of  the  fender  in  question.  The  load  is  con¬ 
sidered  to  be  applied  in  small  increments  associ¬ 
ated  with  the  time  interval.  This  P-D  relation 
is  determined  from  statical  analysis  or  tests.  At 
any  time  t,  the  deflection  D  will  be  calculated 
from  the  displacement  of  the  point  of  contact  C 
and  the  deflection  of  the  structure,  which  will 
equal  to; 


Dx=Xc~X,  ‘ 
Dy  =  Yo—  Yt 


and  D=  •/Dx,+Dyf 

. (7) 


In  considering  the  second  impact,  the  velocity 
of  the  point  of  the  first  contact  at  separation, 
magnitude  and  direction  is  needed,  which  will 
equal  to; 


a,=tan-l  (Xc*!Yc,)t  \ 


(8) 


At  time  of  contact  0^=0,  the  angle  that  the  ap¬ 
proaching  velocity  makes  with  Y-axis,  a  is  con¬ 
sidered  positive  when  the  velocity  vector  of  Y* 
at  G  points  towards  the  point  of  contact  C. 

If  the  value  of  the  spring  constant  of  the  ship 
hull  at  point  of  contact  is  available,  the  elastic 
deformation  of  the  ship’s  hull  can  be  evaluated. 
Let  Kk.  Sk.  and  /u  define  the  spring  constant, 
deflection,  and  the  damping  coefficient  of  the 
hull  at  point  of  contact  respectively,  then  the 
equation  of  deflection  of  the  ship's  hull,  in  P- 
direction,  will  be; 

:  Sk-{P-Xk-Sk-M-Sk)  . (9) 

In  this  case  P  will  be  function  of  (Ac,  Yc.  X$.  Yt, 
Skh •  The  initial  conditions  of  motion,  at  time  of 
contact,  are; 


Xo=Yg=e=p=Xc=Y0 

=AT.=  Y,=S»= 0.0 
Xc~Xc—  V*  sin  a 

Yg—Yc—  V»cos  o 
^  —  p=Xt—Y,=Sk—0.0 


(10) 


The  solution  of  dynamic  equations  is  carried  out 
by  numerical  integration  methods  with  the  help 
of  the  digital  computer. 


(3)  Energy  Equations 

The  developed  previous  equations  are  valid  as 
long  as  the  ship  being  in  contact  with  the  fender. 
During  this  time,  the  following  energy  equations 
are  valid; 
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I.  Part  of  the  ship's  kinetic  energy  trans- 


mitted  to; 

i) 

Fender 

Vf^P-dDi 

ii) 

Structure 

V.^K.-dSt 

iii) 

Ship’s  Hull 

V*= j'  Kx-dS*t 

^effective  energy  Et 

II.  Work  done  by  the  ship  in  rotational 
motions; 

i)  Swinging  Wi  =  ^  P-r-dOt 

ii)  Heeling  P-H>d<h 

111.  Part  of  energy  induced  in  the  system 
vibration; 

i)  Ship  &k.-jAf,(Xoi+to>) 

ii)  Structure  E,=^rM\(Xf+Y J) 

2  4 

. (11) 

The  above  equations  should  satisfy  the  con¬ 
servation  of  energy  during  berthing  i.e. 

Et=—MtW—Et+Eu.+E$  . (12) 

( 4  )  Broadside  Berthing 

If  the  motion  of  the  ship  during  the  berthing 
operation  is  mainly  governed  by  tugboats,  as  is 
always  the  case  with  the  large  ships,  the  ship 
can  make  contact  with  the  berthing  structure 
entirely  broadside.  In  this  case  the  ship  will 
undergo  dynamic  motion  which  has  two-degree 
of-freedom;  swaying  and  rolling,  and  terms  con¬ 
taining  8  in  equations  (5)  and  (11)  with  vanish. 

If  the  energy  dissipated  by  rolling  motion  is 
neglected  for  safety,  then  the  fender  system 
will  be  designed  to  absorb  all  of  the  kinetic 
energy  of  the  ship,  which  is  the  cate  when  C= 
1.0  in  equation  (2).  This  mode  of  berthing  is 
considered  ideal  as  the  ship  impact  load  will  be 
uniformly  distributed  on  the  structure1'). 

If  we  assume  that  the  berthing  structure  it 
provided  with  a  fender  which  hat  a  linear  apring 
constant  Kj  the  equations  of  motion  will  be; 

i)  Structure 

y,=[P-KrY,-v?,)IM: 


ii)  Ship  (+ fender) 

Y,=-PIMt 

in  which  P*Kj(Y,-Y,)  and  Y„Y,  are  the  dis- 
placements  of  the  structure  and  the  point  of 
contact  respectively  in  F-direction  Fig-  1. 

If,  at  the  time  of  contact,  the  following  condi¬ 
tions  exist; 

Y,=  Yt=0.0  F,= 0.0  Yt=V,  . (1<1) 

and  assuming  p— 0.0,  the  analytical  solution  of 
these  eq.  (13)  is  given  by; 

y,=A)  sin  #»i/+ Ai  sin  pit 
Yt-MBi  sin  p,t+AiBi  sin  pit 
The  equations  of  energy  become: 

V,^1  K.Y.dY,^  K.Y.Y.dt  '  . 

=  l/2/f»(Ai  sin  p,t + At  sin  pity 
V>=5‘  KAYc-Y,)(Yt-Y,)dt  . (16) 

=V2KT[AxfBi-l)t\npd 
+Ai(Bj— l)sin  pity  t 

In  the  case  of  a  very  rigid  berth,  which  offers 
a  large  resistance,  the  deflection  Y»  will  be  very 
small  and;  consequently,  its  ability  for  energy 
absorption  will  be  very  poor,  and  can,  therefore, 
be  neglected.  In  this  case  all  the  portion  of 
energy  consumed  in  the  swaying  motion  of  the 
ship  should  be  absorbed  by  the  fenders.  , 

Notations 

/)_!=: polar  moment  of  inertia  about 
the  longitudinal  axis  (1—1) 
passing  through  the  C.C. 

7j_i=polar  moment  of  inertia  about 
the  vertical  axis  (2—2)  pass¬ 
ing  through  the  C.G. 

Afi = the  effective  mass  of  the 
structure. 

Mi=the  virtual  mass  of  the  ship 
while  swaying. 

Fsthe  ship  acting  load. 

H,  Hi—iht  vertical  distances  between 
the  C.G.  and  the  point  of  con¬ 
tact  and  the  meta-center  re¬ 
spectively. 

rsthe  distance  from  the  C.G.  to 
the  point  of  contact. 
w=the  ship's  displacement 

weight.  . <17> 

yacthe  angle  that  the  velocity 
vector  makes  with  the  arm 
r  at  time  of  the  first  contact. 
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K„  Restructure  stiffness  in  X  and 
Y  direction. 

r,  Restructure  damping  coefficients 
in  X  and  Y  direction. 

A',.  ,Yr=the  displacement  of  the  struc¬ 
ture  in  X  and  Y  direction. 
Ai^Vdpi&i-Bt) 

Aj=  —  Vtlpi(Bi—Bt) 
fl,=(a+b)l2*s/l(a+b)l2y+l>c 
Bi^bHpt—o) 

Bt-bHpt'—a) 

a=(K,+K,)IMl 

bmKflMi 

c-KjjMx 

3.  SELECTION  OF  THE  PARAMETERS 
FOR  SOLUTIONS  OF  THE  DYNAMIC 
EQUATIONS 

The  dynamic  equations  of  motion  presented  in 
section  2  contain  the  following  parameters; 
for  the  ship,  Vi,  Mu  h~u  h-u  Hu 
for  the  structure,  Mu  K,, 
for  the  fender,  (P-D) 

and  in  practical  design,  the  ship's  displacement 
weight  (IV)  is,  generally,  the  only  value  given 
with  the  information  for  the  marine  structure  to 
be  constructed. 

This  section  includes  equations,  tables,  and 
graphs  to  help  in  selecting  or  computing  the 
unknown  parameters  for  solutions  of  the  dynam¬ 
ic  equations. 

(1)  Approach  Velocity,  V* 

From  eq.  (1)  the  impact  energy  is  proportion¬ 
al  to  the  square  of  the  approach  velocity.  Thus, 
the  energy  level  will  increase  considerably  if  the 
velocity  is  only  slightly  increased.  In  the  selec¬ 
tion  of  this  velocity  for  design,  many  factors 
should  be  considered,  such  as: 

1)  Method  of  docking 

A  ship  approaching  the  berth  under  the  con¬ 
trol  of  tugboats  usually  berths  with  less  velocity 


than  one  approaching  under  its  own  power. 

2)  Berthing  conditions 

In  the  case  where  the  berth  is  exposed  to  wind, 
waves,  currents,  etc.,  it  is  more  difficult  to  con¬ 
trol  the  ship  velocity  than  in  a  sheltered  berth 
and  the  velocity  may  become  large.  However, 
for  extremely  high  wind  velocities  in  the  order 
of  100  to  120  mph  that  occur  during  short  peri¬ 
ods,  it  is  advisable  to  require  ships  to  temporari¬ 
ly  anchor  away  from  the  berth  in  order  to  avoid 
being  damaged*”. 

3)  Ship  size 

Larger  ships  are  always  berthed  with  great 
care  and  with  the  assistance  of  tugboats.  It  is 
generally  assumed  that  the  larger  the  ship,  the 
smaller  will  be  the  velocity  with  which  the  ship 
will  contact  the  fender  system”. 


Fig.  2  Berthing  velocity  normal  to  dock 
vs.  Ship  .displacement  weight  (after 
Lee).  . 


T.  Lee  presented  the  curves  shown  in  Fig.  2 
from  which  the  berthing  velocity  may  be  select¬ 
ed  for  design.  Under  various  conditions  of  berth¬ 
ing,  Vasco  Costa”  recommended  Table  1  as  a 
guide  for  the  selection  of  velocity. 

Before  proceeding  to  the  selection  of  other 
parameters,  some  relations  concerning  the  ship 
characteristics  will  be  discussed.  For  example, 
a  tanker  of  length  L,  draft  d  and  beam  breadth 
B,  Fig.  1,  will  have  the  following  empirical  re¬ 
lations: 

W=cepLBd  1 

cj=0.75  to  0.80  =  0.78  | 


Table  1  Approach  velocity  of  berthing  ships 

(ft/sec) 


Approach 

Displacement  o(  the  ship 

conditions 

Up  to  3,000  ton 

Up  to  10,000  ton 

Over  30,000  ton 

Strong  wind  and  swell 

Difficult 

2.5 

2.0 

1.5 

Strong  wind  snd  swell 

Favourable 

2.0 

1.5 

1.0 

Moderate  wind  and  swell 

Moderate 

1.5 

1.0 

0.1 

Protected 

Difficult 

1.0 

o.t 

O.S 

Protected 

Favourable 

0.1 

o.s 

0.4 

(slier  V.  Cotta) 
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D/d  stotal  depth/drafts  1.33 1,1 
HD- 13.75 
l/rf«l 8.25 

5S=d-0.S2D*0.31d 

radius  of  gyritioD  about  vertical 
axis  through  G=0.2f.,) 
k,  -  radius  of  gyration  about  longi¬ 
tudinal  axis  through  G=«'2f 
*'« 0.37  to  0.47&0.42 


(2)  Derivation  of  Added  Mass  Equations  for 
Berthing  in  Shallow  Water 

1)  Added  mass  in  horizontal  motion,  Mi 

It  is  well  known  that  when  a  ship  moves  from 
deep  to  shallow  water,  as  in  the  case  when  berth¬ 
ing,  the  added  virtual  mass  is  increased  due  to 
the  presence  of  restricting  boundaries.  Koch1” 
investigated  the  effects  of  shallow  water  on  add¬ 
ed  virtual  weight  for  both  vertical  and  horizontal 
vibration.  The  measurements  were  made  for  a 
block  having  a  half-beam  b  and  a  draft  d.  The 
results  of  the  experiments  are  shown  in  Fig.  3. 
To  apply  the  Koch  results  on  hull  shape  sections, 
it  is  suggested  that  the  added  weight  to  any 
section  of  the  ship,  calculated  for  deep  water  for 
this  particular  section  shape,  should  be  increased 
by  the  ratio  of  the  added  weights  in  shallow  and 
deep  water  for  a  rectangular  section  of  the  cor¬ 
rect  beam-draft  ratio  and  depth  of  water. 

Furthermore,  on  the  basis  that  many  modern 
-niters  have  very  similar  hull  shapes,  sections 


Fig.  3  Effect  of  shallow  water  on  added  mass 
in  horizontal  vibration,  Curves  of  Cm 
for  rectangular  section  (after  Koch). 


Fig.  4  Actual  sections  of  ship  (after  Kumai), 


Fig.  S  Values  of  added  mass  coefficient  Cjr 
for  horizontal  motion  in  deep  water. 


for  typical  tankers,  Fig.  4 »>,  are  considered  in 
the  analysis  hereafter. 

With  the  help  of  the  curves  of  Fig.  5  and  these 
sections,  the  distribution  of  the  added  mass  along 
the  ship  length  may  be  obtained  (step  6  in  Table 
2,  and  curve  a  Fig.  7).  Considering  the  Koch 
results  for  rectangular  sections,  and  the  proce¬ 
dure  above,  the  added  mass  of  actual  hull  shape 
sections,  where  the  ratio  of  water  depth  T  to 
draft  d  is  2.0,  can  be  evaluated  (steps  7  to  10, 
Table  2).  The  results  are  shown  by  curve  b  of 
Fig.  7). 

However,  the  water  depth  considered  by  Koch 
was  deeper  than  that  required  for  berthing.  The 
experimental  results  obtained  by  Marwood  and 
Johnson,  Fig.  6  are  of  considerable  help  in  this 
field,,).  This  is  bearing  on  the  fact  that  the 
percentage  increase  in  Cm  in  shallow  water, 
where  77<f=2.0,  than  that  in  deep  water  for  the 
ship  mid-sections  (4-7),  computed  on  the  basis 


Fig.  $  Effect  of  shallow  water  on  added  mass 
—(Marwood  &  Johnson). 
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Table  2  Evaluation  of  added  mass  in  shallow  water 


of  the  Koch  results,  showed  a  close  agreement 
with  Marwood  and  Johnston’s  experimental  re¬ 
sults,  as  indicated  in  Fig.  6. 

Taking  into  account  the  results  shown  in  Fig. 
6,  the  derived  values  of  the  added  mass  coeffi¬ 
cients  in  step  10  are  re-calculated  for  a  water 
depth  and  draft  ratio  of  1.2  (step  11,  Table  2  and 
curve  c.  Fig.  7). 

By  integrating  along  the  length  and  substitut¬ 
ing  for  </=B/2.61  from  Table  2,  the  following 
equation,  for  the  added  mass  in  horizontal  vibra¬ 
tion,  was  derived; 


Fig-  7  Distribution  of  added  mass  along 
length. 


m’tz^pdBL  (19) 

The  virtual  mass  becomes 

Afj=m+m'  (20) 

and  the  virtual  moment  of  inertia  /s— a,  in  yawing 
motion  will  be°>° 

Is-s—Mik*  (21) 

2)  Checking  the  derived  formula 
Taking  into  account  both  model  and  prototype 
experiments,  Vasco  Costa0  presented  the  follow¬ 
ing  equation  for  estimating  the  virtual  mass  of 


a  berthing  ship; 

M,=M+m'mm(l+2dlB)  . (22) 

Shu  T’ien  Li>»  also  presented  the  following  equa- 
ti  an; 

Af, = w(l + *B/(  1 6D))  . (23) 

The  virtual  mass  for  different  tankers,  using 


the  three  equations,  was  calculated,  Table  3.  It 
can  be  seen  that  the  derived  equations  yield  re¬ 
sults  which  are  about  iff  more  than  the  Shu 
T’ien  Li  equation  and  about  iff  less  than  the 
Vasco  Costa  equation. 

3)  Added  mass  of  inertia  in  rolling  motion 

The  added  mass  moment  of  inertia  coefficients 
depend  upon  the  sectional  shapes  and  the  ratio 
of  beam  to  draft;  the  same  parameters  which 
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Table  3 


D/W 

Ditpl.  W. 

L_ 

m 

s 

i 

1 

11 

111 

111 

I 

Tankers 
10  000 

ton 
13  300 

m 

140 

m 

17.2 

m 

7.9 

ton 

2SS00 

ton 
19  000 

ton 
24  700 

97.0  * 

20  000 

2S700 

ITS 

22.4 

s.s 

41000 

39  000 

47  750 

99.9 

30  000 

40  000 

200 

25.1 

10.3 

72  000 

59  600 

66  200 

95.0 

50  000 

S6500 

230 

32 

11.4 

114  000 

103  000 

105700 

92.5 

*5000 

113  000 

260 

31.1 

14.0 

196  000 

173  000 

180  000 

92.0 

100  000 

133  000 

2S5 

41.2 

14.6 

223  000 

206  700 

213  000 

95.0 

l*Vaic  Cotta  Formula  average:  95% 

1l“Shu  T’ien  Li  Formula 
ll!*The  New  Formula 


Fit:.  8  Distribution  of  added  virtual  weight 
along  length  in  torsional  vibration 
(after  Kumai). 


effect  Cjr.  Also,  these  coefficients  depend  on  the 
location  of  the  center  of  rotation.  Model  experi¬ 
ments  were  carried  out  by  Kumai  on  prismatic 
models  having  sections  corresponding  to  those  of 
a  tankerU).  Applying  the  experimental  results 
on  the  actual  sections  of  a  tanker  ship,  Fig.  4, 
Kumai  obtained  the  distribution  of  the  added 
mass  moment  of  inertia  along  a  typical  tanker 
hull  in  the  loaded  condition.  Fig.  8.  By  integrat¬ 
ing  along  the  length,  Kum»i  derived  the  follow¬ 
ing  expression; 

J/i=0.00531(l +0.365<//d/)S<i  <ton-m») 

. (24) 

For  a  fully  loaded  tanker  (d=dj).  the  above  ex¬ 
pression  becomes; 

d/,=0.00725B«L  (ton-m») 

. (25) 

Eq.  (24)  represents  the  added  mass  moment  of 
inertia  in  deep  water,  but,  as  was  discussed  pre¬ 
viously,  berthing  always  takes  place  in  shallow 
waters.  Therefore,  the  effect  of  shallow  water 
on  the  inertial  moment  baa  to  be  considered. 

Matsuura  and  Kawakami»>  performed  numeri¬ 
cal  computations,  applying  tbe  finite  element 
method,  on  the  effect  of  the  restricted  water  on 
the  inertial  coefficients,  Cr.  Rectangular  sections 
having  a  ratio  of  faalf-beam  to  draft  of  1.0  and 
two  locations  of  center  of  rotation,  yi/«f*=0.0  and 
1.5,  were  investigated,  Fig.  9. 


To  apply  the  results  shown  in  Fig.  9  to  the 
case  of  a  fully  loaded  ship,  the  following  correc¬ 
tions  are  necessary;  a)  first,  since  we  are  inter¬ 
ested  in  the  moment  of  inertia  (/i_i)  about  the 
center  of  gravity  of  the  ship  where  yt/d=0.31, 
eq.  (18),  the  effect  of  shallow  water  on  the  cen- 


Fig.  9  Effect  of  shallow  water  on  added 
mass  moment  of  inertia  in  rolling 
(after  Matsuura  &  Kawakami). 


Fig.  10  Variation  of  Cr  with '  respect  to 
center  of  rotation  for  b/dm  1.0  (after 
Matsuura  &  Kawakami). 
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ter  of  rotation  becomes  important.  This  is  ob¬ 
tained  by  comparing  the  inertia!  coefficient:,  Cr> 
at  point  ot  where  T/4=  1.2,  Fig.  9,  with  point  at 
in  Fig.  10.  That  is, 

Cr«)/Cr*,=0.12/Q.l«=1.2 

From  which  the  value  of  Cr  at  y»/4=0.31  and 
774=1.2,  namely  Cre,  will  approximately  be 
equal  to 

Crs= 1.2(0.056) =0.067  . (26) 

b)  The  second  correction  is  obtained  from  the 
consideration  of  three  dimensional  motion.  The 
correction  factors  from  2  to  3  dimensional  mo¬ 
tion  can  be  obtained  from  Fig.  II1”.  If  L/4= 
18.25,  eq.  (18),  the  correction  factor  correspond¬ 
ing  to  pure  rolling,  (»=0),  is  equal  to  0.96. 
Multiplying  eq.  (26)  by  this  value  yields; 

Cru =0.067(0.96) = 0.064  . (27) 


Fig.  11  Taylor  correction  factor  from  2  to 
3  dimensional  motion. 


c)  The  third  correction  is  obtained  from  the  ef¬ 
fect  of  actual  ship  hull  sections.  A  comparison 
was  made  with  sections  4  to  7,  Fig.  8,  having 
Ct  in  deep  water  equal  to  0.06,  and  the  rectan¬ 
gular-shaped  sections  where  Cr  in  shallow  water 
was  derived,  eq.  (27).  From  Table  2  the  area 
coefficient  (a)  of  these  sections  is  equal  to  0.996 
or  about  1.0,  giving  a  very  slight  effect  due  to 
the  round  edges  of  these  particular  sections. 
The  other  sections  have  already  been  considered 
by  the  use  of  tbe  equation  derived  by  Kumai, 


eq.  (24).  Hence,  for  deduction,  the  Cr  values  of 
ship  sections  in  shallow  water  (values  included 
in  Fig.  8)  can  be  multiplied  by  the  ratio  0.064/0.06, 
which  denotes  the  Cr  value  in  shallow  water  as 
compared  to  deep  water  for  sections  4  to  7.  This 
wilt  lead,  finally,  to  multiplying  eq.  (25)  with  the 
above  ratio  for  obtaining  the  added  mass  moment 
of  inertia,  as  follows,  in  shallow  water; 

4/i=0.00774£<I  (ton-m»)  . (28) 

The  polar  moment  of  inertia,  /»,  about  a  longi¬ 
tudinal  axis  passing  through  C  is  equal  to 
/t=m(fc,)> 

and  substituting  from  eq.  (18)  results 
h—0.7&fiLBd(0.42B')> 

Letting  4=21/2.62  from  Table  2,  step  5,  p=1.03 
’ton/m* 

/«=0.05412?U.  (ton  m») 

Jh!(!t) = (0.007740  <L)/(0.0541£'L) = 0.1 43 

. (29) 

This  ratio  is  in  close  agreement  with  the  ratio 
of  0.15  given  by  Prof.  Hayashi”.  The  virtual 
moment  of  inertia  /i-t  becomes 

/i-i=l.l43m(*,)»=0.202B»m  (ton-m>)  j 
or  /i_i=/,+4/,=0.0618fla  (ton-m>)  j 

.  . (30) 

4)  Hi  or  CM,  which  denotes  the  vertical  dis¬ 
tance  between  the  ship's  center  ot  gravity,  C, 
and  its  metacentre,  M,  Fig.  1,  can  be  calculated 
from  the  following  equation1”; 

B=LICi+CM(d/D)C,  (31) 

where  Cj  and  Ct  are  constants  having  the  fol¬ 
lowing  values  for  oil  tankers  (where  KU/D= 0.52); 

•  Cj=12.5  (U-shape)-~13.2  (V-shape) 

C»=5.7 


Fig.  12  Mitsubishi  Shoji  Oil  Berth. 
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5)  Effective  mats,  Mu  of  the  structure 
For  determining  the  value  of  Mi  that  ahouid  he 
substituted  in  the  motion  equations,  the  authors 
investigated  some  existing  berths  in  Kobe  harbor. 

One  is  an  oil  berth  belonging  to  the  Mitsubishi 
Shoji  Co.  which  consists  of  four  dolphins  of  dif¬ 
ferent  sites  and  a  platform,  Fig.  12,  located  back 
from  the  berth  line.  The  berth  is  provided  with 
V-type  rubber  fenders.  The  berth  was  designed 
to  accomodate  tankers  varying  in  aize  (8  000, 
20  000  and  47  000  DWT)  under  the  control  of 
tugboats. 

a)  The  47  000  DWT  berths  against  two  dol¬ 
phins  which  have  .an  effective  weight  equal  to 

cap+effective  piles  weight 

*=2(160+180)*= 680  tans. 

b)  The  20  000  DWT  should  use  two  dolphins 
which  have  an  effective  weight  equal  to 

2(100+ 120) =440  tons. 

c)  The  8  000  DWT  should  berth  against  two 
dolphins  which  have  an  effective  weight  equal  to 

2(60  +  80)=280  tons. 

The  displacement  weights  of  the  above  ships 
are  approximately  (1.3  W)  or  61  000.  26  000  and 
10  400  tons,  respectively.  The  virtual  masses 
(Afi)  will  be  equal  to  104  000,  48  000  and  25  000 
tons,  respectively,  Table  1  and  the  ratio  MJMi 
will  be  equal  to  680/104  000=0.007,  440/48  000= 
0.009,  and  280/25  000=0.011. 

Mi/Mi=0.01  . (32) 

Thus  for  the  first  design  approximation,  we  can 
assume  that  Mi  will  be  as  much  as  of  the 
virtual  mass.  Mu  of  the  approaching  ship.  If, 
at  the  end  of  the  calculations,  the  difference  be¬ 
tween  the  derived  value  of  Mi  and  the  assumed 
value  is  great,  the  calculations  can  be  repeated 
using  the  derived  value  for  Mu 
6)  Structure  and  fender  resistance 
As  discussed  previously  for  the  design  of  marine 
structures,  the  function  of  the  structure  should 
be  known  in  advance  in  addition  to  the  ship  dis¬ 
placement  weight.  In  the  selection  of  the  resis¬ 
tance  (At)  of  the  structure,  the  structure  func¬ 
tions,  whether  rigid  or  flexible,  should  be  con¬ 
sidered.  If  the  berth  carries  heavy  loads  (heavy, 
delicate  equipment  carried  on  the  deck,  cranes, 
power  station,  etc.),  there  is  no  choice;  deflection 
must  be  limited.  The  construction  should  be 
rigid  and  provided  with  elastic' fenders  to  absorb 
the  ship’s  kinetic  energy.  On  the  other  hand, 
if  deflection  it  allowed,  the  berth  can  be  flexible. 

A  pile  dolphin  is  a  good  example  for  this  type. 


4.  SOLUTIONS  OF  THE  DYNAMIC 
EQUATIONS  BY  NUMERICAL  IN¬ 
TEGRATION 

The  differential  equations  of  motion  of  the  ship 
and  the  berthing  structure  are  in  the  form  of: 

Structure  X,=f(X„  P(X„  Xc),  X.) 

Ship 

i)  roll  f=G(t,P(X„Xc),*i) 

ii)  yaw  i=Gi(6,  P(X„Xt).  i) 

iii)  sway  Xc=Gi(Xc.  l,t.P{X„X<), 

(Xe—Xt))  l 

. (33) 

(These  notations  are  explained  in  section  2.) 

In  eq.  (33)  the  first  differential  terms  represent 
the  damping  effect  which  is  generally  neglected. 
Solutions  of  eq.  (33)  are  quite  difficult  to  be  car¬ 
ried  out  analytically,  especially  in  the  case  where 
fenders  (such  as  rubber)  that  have  non-linear 
spring  constants  are  used. 

A  numerical  solution  implies  the  determination 
of  the  displacement  and  velocity  of  a  system  as 
a  function  of  time.  These  displacements  and 
velocities  are  obtained  in  a  step-by-step  integra¬ 
tion  procedure,  starting  with  given  initial  condi¬ 
tions.  There  are  many  different  methods  of 
numerical  integration  from  which  two  methods, 
the  Newmark  p  method’*’  and  the  Runge-Kutta- 
Gill  method’”,  will  be  explained. 

( 1 )  Time  Interval  Effect  on  the  Two  Numeri¬ 
cal  Integration  Methods 

Tests  have  been  conducted  to  study  the  effect 
of  the  time  interval  on  the  accuracy  of  the  two 
methods.  In  these  tests,  fenders  with  linear 
spring  constants  were  used,  for  which  the  exact 
solution  of  the  equation  of  motion  is  obtained 
by  using  eq.  (15).  Through  the  comparison  of 
results  which  are  included  in  Table  4  and  Fig. 
13,  the  following  conclusions  could  be  made: 

1)  The  error  involved  in  fender  absorbed  ener¬ 
gy,  V>,  is  relatively  small  compared  to  the 
energy  V,  absorbed  by  the  structure. 

2)  The  percentage  of  error  involved  in  the  V$ 
values  is  nearly  twice  that  of  the  structure 
maximum  deflection,  Xu  >.c.  dViJlf=2 JXifi. 

3)  For  the  same  time  interval  the  error  in¬ 
curred  using  the  Runge-Kutta-Gill  method 
is  greater  than  that  involved  in  the  New¬ 
mark  p  method.  The  difference  also  in¬ 
creases  as  the  time  interval  is  increased. 
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Fir.  13  Convergence  of  error  with  respect 
to  time  interval. 


(2)  Selection  of  a  Suitable  Time  Interval 
Through  many  investigations  carried  out  by 
the  authors,  Fig.«4,  a  formula  linking  the  time 


Fig.  14  Effect  of  time  interval  on  the  berth 
of  berthing. 


interval  variation,  the  berthing  data,  and  the 
error  involved  in  the  maximum  deflection  of  the 
berthing  structure  was  developed. 

The  curves  shown  in  Fig.  15  were  plotted  from 
Fig.  14  for  an  error  equal  to  2#  of  the  struc¬ 
ture  maximum  deflection.  This  gives  2£  error 
in  the  structure  maximum  reaction  (KtX ,)  and 
nearly  Afi  in  the  structure's  stored  energy.  In 
this  figure  the  rt-axis  represents  the  berthing 
data,  which  is  the  factor  N,  and  the  y-axis  re. 
presents  the  time  interval  At.  The  factor  N  is 
a  function  of  the  ship's  virtual  mass,  the  struc¬ 
ture’s  effective  mass,  the  structure  and  fender 
spring  constants,  and  the  ship's  approaching 
velocity  according  to  the  following  equation; 

A^fM/AMiXlOO)  (cm)  . (34) 

where  At  is  obtained  from  eq-.  (17). 

From  Fig.  15  the  relation  between  N  and  the 
time  interval  At  for  an  error  of  2V  in  AXi  is 
given  by: 

i)  The  Newmark  0  method  (0=1/4) 

At = 0.00067 N + 0.0012  (sec).  . (35) 

A' £3.73  ■ 

df=0.00015AM— 0.0029A* 

^  ■ 

+0.018A!  —0.0313  (sec)  --(36) 
Af£3.73 

ii)  The  Runge-Kutta-Gill  method 
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its 


Jt= 0.0018  (sec) '  " . (37) 

#£3.73 

d/=0.000066#»-0.00137#»  ' 

+0.0095# -0.0178  (tec)  . (38) 

#2:3.73  -  . 

( 3  )  Application  to  the  Nonlinear  Spring  Con- 
ttar.t  Fen  den 

Both  rubber  and  retractable  fenders  possets 
non-linear  relations  between  the  load  and  the 
displacement.  To  apply  the  preceding  equations 
for  selecting  the  time  interval  for  structures 
provided  with  these  type  of  fenders,  the  load  is 
considered  to  be  applied  in  small  increments  as¬ 
sociated  with  the  time  interval.  The  procedure 
for  calculation  is  as  follows:  ,  .. 

1)  The  spring  constant  [Ks)  of  the  fender  corre¬ 
sponding  to  aero  displacement  is  taken  from 
the  given  load-displacement  relationship. 

2)  The  factor  #  can  be  calculated  and,  con¬ 
sequently,  the  time  interval  from  eqs.  (35), 
(36)  or  eqs.  (37),  (38).' 

3)  Substituting  in  the  equations  of  motion, 
the  displacement  of  the  fender  at  the  end 
of  the  interval  can  be  obtained. 

4)  The  spring  constant  corresponding  to  this 
displacement, can  be  calculated  as  indicated 
in  step  1). 

5)  Repeating  steps  2),  3),  and  4)'  until  the  al¬ 
lowable  fender  displacement  is  reached. 


5.  APPLICATION  "  ■  '  ' 

( 1 )  A  Case  of  Genera)  Berthing  and  Fenders 
of  Non-linear  Spring  Constant 
Soloed  Example 

a)  Data  Given; 

A  tanker  ship  of  displacement  weight  W 
=40  000  ton 

Approaching  velocity  V»= 15  em/sec 
Ship  characteristics  are; 

Length=200  m  Breadth=25.8m 
Draft = 10.8  m 
Berthing  data; 

r=80m  H- 3.0  m  y=50*  «=20* 

Assuming,  for  the  first  trial,  the  fender  system 
data  as; 

A*=15  ton/cm  Fender=2  pieces  of  type  1, 
Fig.  17. 

b)  By  applying  the  formulas  included  in  sec¬ 
tion  3,  the  following  data  was  computed; 

Afi*=70.69  ton-sec*/cm 


Fig.  16  Flow  chart  for  design  procedure  by 
using  newmark  "S"  method. 
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Table  5 


Case 

No. 

Mode  of  berthins 

r* 

r 

m 

H 

m  : 

K, 

ton/cm 

Fender 
Fis.  17 

Effective 

energy 

E, 

c~§«. 

Ei 

% 

# 

i 

tee** 

1 

General  berthins 

so 

so 

2 

IS 

•* 

Type  1 

41.79 

52.4 

0.1112 

xio-* 

0.151 

X10-* 

2 

No-rolling 

# 

m 

0 

IS 

# 

42.09 

53.44 

\ ! 

^  ••• 

so 

to 

0 

25 

4.25 

34.70 

54.5 

3 

upring  constant 

A>*  500  ton/m 
mCr 

37.40 

57.7 

Tablet  3  and  4 
refer.  4) 

*•  The  fender  at  iff  nett  uled  in  these  catet  it  chosen  by  trUlt  to  give  T',"0.0  at  max.  nway.' 
***  A/i“120ton-eee*/em,  10 cm/tec  . 

•“*  Xr“the  comboond  atlffnett  of  the  fender  *yttem»/f»A//CX«+A/)  ■  — *’■  •  - 


Mi— 0.7  ton-sec’/cm 
Hi =577=2-27  m 
/i_i=54  x  10<  ton-sec*-em* 

/*_i=113x  10’  ton-see*-cm* 

c)  For  numerical  integration,  Newmark  p 
method  with  ^=1/4  was  used.  Besides,  equations 
(35)  and  (36)  were  applied  for  selecting  the  time 
interval.  Calculations  were  carried  out  by  the 
digital  computer,  the  flow  of  computations  is 
shown  by  the  block  diagram  Fig.  16.  Results  is 
included  in  Table  5. 

(2)  Besides,  for  comparing  the  presented 
method  with  other  investigators  methods,  two 
other  examples  were  tried.  Table  5. 

6.  COMMENTS  ON  THE  RESULTS 

(1)  Verification  of  the  Developed  Method  of 

.  Analysis  • 

To  verify  the  assumptions  presented  in  estab¬ 
lishing  the  dynamic  equations  included  in  section 
2,  conservation  of  energy  before  and  after  colli¬ 
sion  is  checked  from  the  deduced  results. 

1)  Conservation  of  the  kinetic  energy 

From  Table  5  we  have  at  t^t=0.09t=0.0  . 

i)  Part  of  the  ship's  kinetic  energy  transmit¬ 
ted  to: '  •  ■  *■  '  *" 

Fender  system™ Effective  energy 

£,=41.79  (ton /ro) 

ii)  Part  of  energy  induced  by  the  system  vi¬ 
bration:  ■  —  •• 

Ship  - - 

'  i 

•  £»*jM,<Xo>+i'oV 

+y/i-«^*+,|-/j-i^* 

at  Vo* 0.0  this  eq.  yields  to; 


substituting  'for  6  and  ft  their  values, 
0.11124x10-*  and  0.151  x  10-*  sec-*  re- 
spectively,  we  obtain;'* 

Eu, ** 29.96+6.69+0.616= 37-26  (ton-m) 
Structure  ..... 

£»=-£- Afi-S*  ,.i 

■  1 

x  0.7  x  (— 4.659)*/100=0.077 

(ton-m) 

'  £,*+£.=37.26+0.077=37.3 '  (ton-m) 

Total  energy  at  Vc=0.0  will  be"}' 

ii+£«s+£’.=79.09  (ton-m)  . . (39) 

But  the  ship's  approaching  energy  is; 

/  £•= j x 70.69 x  155/100=  79.520  (ton-m) 

'  -  -t . (40) 

Eq.  (39)=Eq.  (40)  which  satisfied  the  conserva¬ 
tion  of  kinetic  energy  before  and  collision. 

2)  Conservation  of  the  moment  of  momentum 
i)  At  time  when  the  ship  made  the  first  con¬ 
tact,  the  moment  of  momentum  is  given 
by;  '  '  •  •  7-  •'  • 

*70.69'  x  15  x  8  OOOx  0.766/100 
=649.8x10*  (ton-sec-m). 

"■*;>  :  ‘ . (41) 

ii)  At  time  when  the  ship  started  to  rebound, 
the  moment  of  momentum  is  equal  to; 

.  MMi*Iist+Ii-ii+Mi*JXoi+'$'e>  *r 
«126.10xl0*+8.15xl0* 
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+527.80x101 
=662.05x10*  (ton-sec-m) 

. (42) 


7.  COMPARING  THE  METHOD  WITH 
OTHER  INVESTIGATORS  METHOD 


Eq.  (41)=Eq.  (42)  which  satisfied  the  moment 
of  momentum  principle*. 


( 2  )  Checkin*  of  Time  Interval  Selection  Proce¬ 
dure 

The  accuracy  of  the  deduced  results,  as  teen 
in  the  previous  paragraph,  verified,  on  one  hand, 
the  right  procedure  of  developing  the  equations 
of  motion,  and  on  the  other  hand,  it  supported 
the  author's  recommendations  for  the  evaluation 
of  the  time  interval  dT  presented  in  section  4. 
Pig.  18  shows  the  variation  of  dT  through  the 
numerical  integration  process  acoording  to  the 
variation  of  the  fender's  stiffness  shown  by  Pig. 
17.  Any  mis-choice  of  the  time  interval  will  lead 
to  large  errors,  and  sometime,  lead*  to  unreason, 
able  results  as  included  in  Table  4. 


Fig.  17  The  load  and  energy  vs.  deflection 
for  V600H  rubber  fender  (Tokyo 
Rubber  Dock  Fenders). 


Fig.  18  Variation  of  time  interval  dT  vs. 
time  of  berthing. 


i)  Case  2,  Table  5,  in  which  N= 0.0,  is  simi¬ 
lar  to  the  case  treated  by  Vasco  Costa,  i.e.  sway¬ 
ing  and  yawing  motion*  are  only  considered. 
From  the  table  we  have;  at  Vc=0.69+0.0  cm/sec. 

£.=(!'«+  V*+  V/)=(18.40+3. 42+20.27) 

=42.09  (ton/m)  . (43) 

Substituting  with  the  given  data  in  the  equation 
given  by  Vasco  Costa,  we  obtain; 

=42.15  (ton/m)  ' . -(44) 

ii)  Case  3  is  similar  to  the  case  treated  by 
Hayashi  &  Shirai,  i.e.  fender  of  linear  spring 
constant  and  y,  Fig.  1,  equals  to  90*. 

From  i)  and  ii),  if  we  considered  some  error 
due  to  numerical  integration,  the  author's  meth¬ 
od  will  be  in  agreement  with  the  two  apedal 
cases  treated  in  references  4)  and  5). 

8.  SUMMARY  AND  CONCLUSIONS 

The  foregoing  study  describes  an  analytical  - 
treatment  of  the  ship  berthing  problem,  based 
on  the  dynamic  response  of  the  ship  and  the 
fennder  system  during  berthing.  1 

The  presented  analysis  covers  almost  the  main 
factors  involving  in  berthing  operations.  These 
are  comprised  of; 

i)  The  approaching  mode  of  the  ship  with 
reference  to  the  face  of  the  berth,  de- 

.  signed  by  a  and  y. 

ii)  The  location  of  the  point  of  contact  on 
the  ship's  hull,  denoted  by  r  and  H. 

iii)  The  structure  stiffness  and  the  fender 
stiffness,  whether  the  latter  of  a  linear  or 
non-linear  spring  constant,  beside  to  the 
hull  stiffness  at  the  point  of  contact,  if 
available. 

iv)  The  mechanical  behaviour  of  the  fender 
system  with  respect  to  the  variation  of 
the  acting  load  direction.  This  is  very 
important  in  case  of  rubber-like  fenders, 
as  their  energy  absorbing  capacity  is  a 
function  of  the  load  direction. 

v)  Consideration  of  shallow  water  effect  in 
swaying,  yawing  and  rolling. 

For  solving  the  developed  equation*  of  motions, 
recommendations  and  formulas  for  estimation 
and  selection  of  the  different  parameters,  parti¬ 
cularly  the  time  interval,  are  presented. 
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The  solution  of  these  equation*  will  lead  to 

the  evaluation  of  the  data  required  for  design¬ 
ing  the  fender  system; 

i)  Evaluation  of  the  energy  absorbed  by  the 
fender. 

ii)  Evaluation  of  the  energy  absorbed  by  the 
fender  structure  and  its  dynamic  reaction, 
magnitude  and  direction. 
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